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Abstract 



An attempt is made to interpret the various existing experimental 
data on the single spin asymmetries in inclusive pion production by 
""q^' the polarized proton and antiproton beams. As the basis of analysis 

^ , the chromo-magnetic string model is used. A whole measured kine- 

Qj \ matic region is covered. The successes and fails of such approach are 

• outlined. The possible improvements of model are discussed. 



Introduction 

The growing amount of the experimental data on the single spin asymme- 
try (SSA) attracts attentions of the high energy spin community. This is 
understandable since the perturbative quantum chromo-dynamics (pQCD) 
predicts a zero spin effect at the asymptotic energy, while the experiments 
furnish up to now the opposite signatures. For example, SSA shows to be 
nonzero at the RHIC energy of a/s = 200 GeV p. Therefore the different 
approaches, not so strongly justified as pQCD, but able to explain all SSA 
effects are needed. 

In 1988 the model was proposed [2| for interpretation of the SSA in inclusive 
pion productions and A inclusive polarization. The model was able to de- 
scribe the scarce data existing at that time. Those experimental data came 
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up from the unpolarized beams striking the polarized target and such exper- 
iments suffer from the poor statistical accuracies. Since later a lot of new 
data were published and mostly from the interaction of the polarized beams 
with pure hydrogen target. Recently the polarized Relativistic Heavy Ion 
Collider (pRHIC) becomes operational and the first data on the SSA for in- 
clusive were published jT]. The PROZA experiment at IHEP, Protvino, 
produced for the first time the data on SSA in the polarized proton target 
fragmentation region using unpolarized proton and pion beams. 
The main goal of this article is to apply the chromo-magnetic string model 
(CMSM) of the paper [2] to the new experimental data on the inclusive pion 
asymmetries. The CMSM offers a simple analytical expression for SSA de- 
scription. The original paper was devoted mostly to the interpretation of the 
SSA in inclusive 7r° production. In current paper we extend the CMSM to 
the description of the SSA in the inclusive charged pion production too. 
In section 1 we remind the main features of CMSM according to The for- 
mulae of CMSM are specified for two distinct kinematic regions: the central 
(CR) and beam fragmentation (BF) region. Section 2 describes the applica- 
tion of the model to the CR, while the section 3 is devoted to the BF region. 
Also in Section 3 we expand the model to the SSA in the inclusive charged 
pion productions. The unique data of E704 experiment on the SSA obtained 
with the polarized antiproton beam are analyzed too in scope of the isotopic 
spin and charge conservations. In summary we outline the successes and fails 
of the CMSM in description of the SSA in inclusive pion production at high 
energy. We make some predictions and indicate the possible improvements 
of the CMSM model. 



1 The chromo-magnetic string model 

The basis of this model is following. After collision of two hadrons the color 
tube (string) is stretched between them. In simple case this tube contains 
the flux of the chromo-electric field. But such a system is not stable. In 
order to make it stable the chromo-magnetic field should circulate around 
this tube. The interaction of this field with the color magnetic moment of 
the polarized quark leads to a kick transverse to the string axis. This kick 
depends on the spin orientation of the polarized quark. The estimate of this 
kick gives a magnitude of order 6pT — 0.1 GeV/c j2]. If we introduce the 
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invariant differential cross section p = E-^ for polarized quark asymmetry 
we get the relation 

where (+) and (-) arguments in p mean the up and down direction of the 
quark polarization. Since 6pT < Pr, where pt is a transverse momentum of 
the final pion to which the polarized quark fragments we get the following 
left-right asymmetry in the quark emission 

^ , . _ da{pT + 5pt) - da{pT - Spr) _ 
^ da{pT + Spx) + da{pT — Spx) 

opT d-^p d-^p 

Here B is a slope parameter of p as it is defined in standard way. 

6 da da ^ /, n 

The general formula for description of the SSA can be expressed in the fol- 
lowing form 

An{x) = Pg{x)- Ag{x)-w{x). (4) 

Here Pq{x) is a polarization of an initial quark, carrying a portion x of the 
initial polarized proton momentum, Aq{x) is so called a quark "analyzing 
power" defined above w(x) presents a fractional contribution of a chan- 
nel of interest in parton-parton interactions. 

The transverse quark polarization (transversity) , Pq{x), is related to the spin 
dependent structure function g2{x) measured in deep inelastic scattering 
(DIS). Unfortunately this function has no parton interpretation. For the 
100% polarized on mass-shell quark g2{x) = 0. On the other hand it looks 
reasonable to expect that the x-behavior of Pq{x) is similar to that of g2{x)- 
Note that the present data |3| indicate the threshold like form of the func- 
tion g2{x) which becomes very small for x < 0.2. Similar effect is seen in 
experimental data for SSA in beam fragmentation region (see next chapters). 
Therefore such effect should be included in a more precise parametrization 
of Pq{x) in future. Since at moment the transversity can not be directly ex- 
tracted from experimental data one needs the theoretical inputs. According 
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to the non-relativistic quark model (NRQM) the transversity may be taken 
in the following way 

2 1 
Pq{x) = - ■ a; for n^, Pq{x) = - ■ x for 7r° and 

Pg(x) = -i ■ a; for TT". (5) 

So the NRQM predicts distinct effects for the pion SSA. First, the sign of 
asymmetry is positive for vr^ and for 7r°, while it is negative for vr". Second, 
the magnitude of asymmetry is expected to be higher for the vr"*", than for 
7r° and 7r~. These predictions are in general compatible with experimental 
data, though some deviations may happen. 

Another theoretical inputs can be taken from paper |1] which aimed to ex- 
plain the " spin crisis" by accounting the pion cloud around the proton. Start- 
ing with NRQM, after account for pion loops correction one gets A'u=0.79 
(the sum for two u-quarks), and A(i=-0.31 corresponding to Pd = —O.SPu- 
In such a case we can take Pu{x) = 0.7x and Pd{x) = — 0.55x leading to 

Pg{x) = 0.7 ■ X for 7r+, Pq{x) = 0.28 ■ x for 7r° and 

Pg(x) = -0.55-a; for TT". (6) 

Comparing © to (jSJ one can see that the difference is essential only for vr" 
and this fact will be taken into account at comparisons of the calculated SSA 
with experimental data. 

An{x) may depend on Xi;' which is a Feinman parameter, on pt which is a 
transverse momentum of pion, and on s which is a square of the total energy 
in the center of mass system of colliding particles. 

The invariant differential cross section is assumed to be measured before or 
simultaneously with SSA. 

The general expression for weighting factor depends on the parton distri- 
bution function V^°(x) and parton fragmentation function D°'{x) (a-parton 
flavor). Assuming that the polarization information may be carried mostly 
through the polarized quark q and much less by gluons we get the weighting 
(or asymmetry diluting) factor through the following expression 

Assuming that the main contribution comes from qg' (for cr(g)) and from gg' 
(for (y{g)) scattering via the t-channel gluon exchange between the gluon g' 
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from unpolarized nucleoli and quark q or gluon g from the polarized proton 
the contributions from quark and gluon may be written as 

a{q)^Cp /V^(-)D''(z)-, 

J X Z Z 

a{g)^CA V^{-)D^{z)-, (8) 

J X Z Z 

where the color factors = | and Ca = 3, z is a fraction of the polarized 
quark momentum carried by the final hadron. 

The parton and gluon distribution functions were taken in the forms 



V'i{x) =x-v{x) = 2.8-v^-(l 



X 



i2 



V^{x) = x- g{x) = ?,.^-{\-xf. (9) 

The fragmentation functions for three types of inclusively produced pions 
were taken in the forms: 

DlAz) = ^(1 - D-Az) = ^(1 - z), DUz) ^ 0; 

Dl4z) = Dl4z) = Dlo{z) = {l~zy. (10) 

The conservation of the isotopic spin and charge in parton interaction leads 
to the relations 

D:,iz) = DU^)=Dt-i^) = DU^), 

Dl.{z) = Dl,{z) = Di,{z) = DUz)- (11) 

The practical applications of these relations will be given below in appropriate 
places. 

The weighting factor ((7j) was calculated using the relations ((HI), ©, and (fTUj) . 
The final version of function w(x) for different pions may be approximated 
by 



w(x) = — , (12) 

where c is: c+ = 0.48 for tt'^, cq = 0.64 for tt^, and c_ = 0.96 for tt". 

As it is seen from relations (0} and ()12|) the analyzing power is calculated 

from p and w{x) without any free parameter. 

At lower energy, small pt and xp one may expect the significant contributions 
to the pion yields from resonance productions Therefore we arbitrary 
selected experimental data with ^/s > 10 GeV to avoid in some extent this 
problem. 
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2 Central region 

At present we know 3 pieces of the published data on the SSA for inclusive 
tt" production in the central region (no data for charged pions). They are: 

1. Inclusive SSA for reaction 

p(|)+p^^o + X . (13) 

Experiment E704 at Fermilab measured this reaction at 200 GeV/c. 
Kinematic region: -0.15 < xp < 0.15, 1.48 < priGeV/c) < 4.31 0. 

2. Inclusive SSA for reaction 

p(t)+p^7r° + X . (14) 

Experiment E704 at Fermilab measured this reaction at 200 GeV/c. 
Kinematic region: -0.15 < < 0.15, 1.48 < priGeV/c) < 3.35 

3. Inclusive SSA for reaction 

p + p(T) ^7r° + X . (15) 

Experiment PROZA-M at IHEP, Protvino, measured this reaction at 70 
GeV/c. Kinematic region: —0.15 < xp < 0.15, 1.05 < pT{GeV/c) < 

2.74 m. 

In following we discuss these data in listed above order. 
1. In E704 experiment up or down (to the horizontal plane) polarized proton 
beam of 200 GeV/c striked the unpolarized proton target. The inclusively 
produced vr'^-s were detected by the electromagnetic calorimeters around the 
central region, that is, at the emission angle of ~ 90° in c.m.s. Since the 
experimental invariant differential cross section, p, is almost exponential in 
function of pt its slope parameter B is practically constant. For example, 
at 200 GeV/c the slope parameter B=(4.19±0.08)(Ge1//c)-^ [6J. Then the 
quark analyzing power, Aq = O.IB = 0.419. 

According to the previous discussions © and Q we take the transversity 
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Pq{x) = |x, constant factor co=0.64 p2|l and inserting all above relations 
into we got the final result for 7r° asymmetry in the central region 



0.14 -x^-^ 



(16) 



v/x + 0.64(1 -x) 



4.5 ■ 



One can see several consequences of this formula ()16|) . First, it is scaled 
on argument x= xt = ^ assuming that the slope parameter is constant. 
Second, it becomes zero at x=0 mostly due to a quark polarization and in 
some extent due to the gluon contribution (second term in denominator). 
Third, when x approaches to one asymmetry increases as linear function of 
X. All these behaviors remind the corresponding experimental asymmetry. 
The result of the model prediction according to the formula (|TBjl is presented 
in Fig.l by the solid line. This prediction is consistent with the experimental 
data at 200 GeV/c. 

2. The second group of data is relevant to the SSA in inclusive 7r° production 
by the polarized antiproton beam of 200 GeV/c (E704 experiment) (|T^ [H]. 
This reaction differs from (fT^ by charge of quarks, that is, all quarks in 
polarized proton are replaced by the corresponding antiquarks in polarized 
antiproton. Therefore the charge conservation requires the equality of the 
SSA's in these two reactions, namely 



As seen (Fig. 1) this conclusion is consistent with the experimental data in 
frame of error bars. 

3. Third piece of new data |7| is relevant to the reaction (fT3j) . The data were 
taken at 70 GeV/c unpolarized proton beam with polarized proton target. 
They are presented in Fig.l with opposite sign of asymmetry (to make a 
comparison to the polarized beam case). The expression for asymmetry at 
70 GeV/c deduced from the CMSM looks like 



In comparison to (fTBj) there are two changes. First is relevant to the energy 
difference and second- to the different slope parameters. At 70 GeV/c this 
parameter B=(5.89 ± 0.08)G'eV^~^ according to [Zj. The result of calculation 
is presented in Fig.l by the dashed line. In this case one can see a good 



ANm) +P^n' + X) = A^(p(T) +p^n^ + X). 



(17) 




(18) 
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Figure 1: Asymmetry in function of "Pt in the central region for reaction 
p + p ^ 7r° + X. At 70 GeV/c the target is polarized, while at 200 GeV/c 
the beam is polarized. SSA is also shown for reaction + p ^ 7r° + X 
at 200 GeV/c. The model predictions are ploted by the sohd hne for 200 
GeV/c and by the dashed hne for 70 GeV/c. 
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consistency of prediction and experimental data. In order to check the quan- 
titative predictions of the CMSM one needs much better statistics. 
4. p(t)+P ^ 7r°+X, = 200 GeV (PHENIX experiment at RHIC). At mo- 
ment there are no data on the analyzing power of this reaction at ^/s = 200 
GeV. But PHENIX collaboration recently published the precise data on the 
invariant differential cross section of this reaction for 1 < ^ 13 GeV/c 
region at pseudorapidity | |< 0.39 |HI. This cross section was parameterized 
by authors in the following form 

Here A = 386 mb ■ GeV'^, po=l-219 GeV/c and n=9.99. Now we can 
calculate the slope parameter for the PHENIX invariant differential cross 
section 

BM = " (20) 

dpT Po + Pt 

It is seen from (j20j) that B{pt) decreases with growth of px as p^^ and leads 
to the decrease of asymmetry at large transverse momentum. This trend is 
similar to one stemming out at the lower initial momentum 200 GeV/c) 
in]. The p at momenta 100 and 200 GeV/c was presented in the following 
form 

p = E^ocip'x + MY.{l-XTf. (21) 

The fit to the experimental data leads to the values: = —5.4 ± 0.2, = 
(2.3 ± 0.3) GeV^ and F = 7.1 ± 0.4, xt = Therefore one can find the 
slope parameter 

w ^ 2A^Pt 2F , , 

"'f"-' = WTW - TiO^- '''' 



The slope parameter at PHENIX ()20|1 decreases in the same way in function 
of Pt as the slope parameter at around E704 energy (j^^ indicating on the 
independence of the slope parameter from the initial energy at -\/s=20-200 
GeV energy region. 

One can make the prediction for asymmetry in the above reaction for PHENIX. 
It seems very small, less than 1%. We hope to see soon the experimental data 
from RHIC on A;v(vr°). 

Concluding the discussions of the SSA in central region we note that better 
precisions in of order < 1% are needed for testing the CMSM. 
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3 Beam fragmentation region 



In this domain at energies ^/s > 10 GeV there are the following pieces of 
data: 

1. The STAR result on the SSA in reaction p{]) n^ + X at ^/s=20Q 
GeV/c. Kinematic domain of measurement: 0.18 < a;^ < 0.59, 1.5 < 
PriGeV/c) < 2.3, | r] \^ 3.8 p. 

2. The E704 resuh on the SSA in reaction p(t) + p ^ 7r° + X at 200 
GeV/c. Kinematic domain: 0.03 < < 0.9, 0.5 < priGeV/c) < 2.0 

m- 

3. The E704 resuh on the SSA in reaction p(t) + p ^ 7r° + X at 200 
GeV/c. Kinematic domain: 0.03 < a;^ < 0.67, 0.5 < priGeV/c) < 2.0 

ini- 

4. The E704 results on the SSA in reaction p{]) + p ^ 7r=^ + X at 200 
GeV/c . Kinematic domain of measurements: 0.2 < xp < 0.9, 0.2 < 
PT{GeV/c) < 1.5|I3. 

5. The E704 results on the SSA in reaction p | +p — > vr^ + X at 200 
GeV/c. Kinematic domain of measurements: 0.2 < xp < 0.9, 0.2 < 
PriGeV/c) < 1.5 13;. 

All listed above experimental data are shown in Fig. 2 and Fig.3. 
The expression for asymmetry Q is still correct, as well as for the weighting 
factor w(x). The main change consists in replacing the argument x by Xi;', 
which presents now the portion of momentum acquired by hadron in the 
fragmentation of polarized quark. SSA may be written in the following way 
(see ©) 

ANixp,pT,s) = Pq{xp) ■ Ag{pT,xp,s) ■ w{xp,c). (23) 
The factor c was defined earlier fsee p2|) ). 

In order to calculate Aq we must know the slope parameter B defined earlier 
(jni) for each reaction at the appropriate kinematic domain. In following we 
discuss point by point the above experimental data on SSA's in the beam 
fragmentation region. 

1. STAR made a measurement of SSA and p (invariant differential cross- 
section) at 1/5=200 GeV and the pseudorapidity | |= 3.8. In order to 
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extract the slope parameter the p measured at the same experiment as func- 
tion oi xp was presented as the exponential function of pt and a fit to the 
data resulted in the slope parameter B=(5.94 ± 0.016) GeV~^ with = 27 
at number of points =5. Sure the fit may be improved by adding a new 
parameter quadratic in pt but for our goal the linear approximation in pt is 
acceptable. 

So the final expression for vr'^ SSA at y^=200 GeV looks like 



The result of calculation is shown in Fig. 2 by the dotted line and it is con- 
sistent with the first STAR experimental data. 

2. Since the E704 collaboration did not make the direct measurement of the 
invariant differential cross section simultaneously with SSA measurement at 
200 GeV/c we took such data from paper In this paper the p was mea- 
sured in the interval of 50-400 GeV/c initial proton momentum and for the 
different laboratory angle of the photon emission 30 < 6L{mrad) < 275. The 
measured photon spectrum was converted to p(7r°) by applying the Stern- 
heimer's method. The transverse momenta of 7r° were in interval 0.3-4.0 
GeV/c. The closest to our kinematic domain are data at laboratory 7r° pro- 
duction angle ^i,=30 mrad. We took data in the region 0.6 < pT{GeV/c) < 2. 
By using the exponential function with slope parameter B we found by fit 
that B = (4.75 ± 0.07) (GeV/c)-^. For 9 experimental points xVd-o.f.=0.7. 
Using the relation xp = 2pT/{\fs-tan9cms) the xp region for experiment was 
defined as 0.03-0.27. In frame of these data it is impossible to go farther in 

Xp- 

Taking the quark polarization as Pq{xp) = (1/3) ■ xp the final formula for 
SSA at 200 GeV/c (E704 experiment) for 7r° is 

^n{xp) = -jr, ^-^^ ^- (25) 

The result of calculation is presented in Fig. 2 as a solid line. One can see a 
quahtative consistency between the model prediction and the data of E704 
experiment. 

3. In the same Fig. 2 we present the SSA data for reaction p(|) +p tt^ + X 
at 200 GeV/c. The model prediction is simplified by using the relation (jl7|) 
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Figure 2: Asymmetry in function of in the polarized beam fragmentation 
region for reaction + p — 7r° + X at \/i=19.4 and 200 GeV. The data 
for reaction + p — > tt'' + X are also shown. The solid line is the CMSM 
prediction for y^=19.4 GeV data, the dotted line is for y^=200 GeV data. 
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which is a consequence of transformation from quarks to antiquarks. As seen 
from Fig. 2 there is a good consistency of model prediction (sohd hne) with 
experimental data. 

4. Figure 3 presents the experimental data from E704 experiment on the SSA 
for the inclusively produced pions at 200 GeV/c in the following reactions: 

P(T) + X, p(t) +p-^ 7r±'° + X. (26) 

For calculation of the quark analyzing power, Aq, we must estimate the 
slope parameters B"*" for positive and B~ for negative pion productions. This 
task was done by using the data from paper pLSj. The invariant differential 
cross-section was taken at a/s = 31 GeV, closest to E704 data. The slope 
parameters were estimated by taking two values of pt =0.82 GeV/c and pt 
= 1.42 GeV/c and averaged over region 0.35 < xp < 0.65. Finally we took 
-8^=5.55 GeV~^ and i?~=5.33 GeV^^. So the expression for 7?+ SSA taking 
Pq{x) = (2/3) ■ X looks like 



^^"^"'^^ " x%' + 0.48(?- x^)4-5- ^2^) 



Prediction of this formula for SSA in inclusive vr^ production by polarized 
protons is presented in Fig. 3a by the solid line. Data (closed circles) confirm 
this prediction fairly well excluding the region xp <0.4. 
Fig. 3b shows the SSA for inclusive vr" production in pp collisions at 200 
GeV/c. According to the CMSM taking Pq{x) = —(1/3) ■ x the following 
relation for tt^ SSA can be written 

A'^~ixp) = - os^n'orr? ^4 5 " (28) 
Xp + 0.96(1 — xp) 

Following the non-relativistic quark model we assumed the negative sign of 
the d-quark polarization. As seen from Fig. 3b the CMSM prediction (solid 
line) follows qualitatively the experimental data (closed circles) but it does 
not give a quantitative agreement. A better description of SSA at the large 
Xp furnishes the dotted line corresponding to the relation Pq{x) = —0.55a; 
stemming out from paper |4j. But in both cases the model needs some 
improvements in order to describe better the smallest xp region. 
5. Fig. 3c presents the SSA versus xp for inclusive 7r° productions by the 
polarized proton (closed circles) and antiproton (open circles) beams. As 
seen from this Fig. 3c the charge conjugation rule (jl7j) is well respected by the 



13 



experiment. For SSA in the inclusive charged pion productions the following 
charge conjugation rule should be respected 

ANiPit) +p^n^+X)^ AMt) + P vr^ + X). (29) 

According to this relation one can predict the inclusive charged pion 
asymmetries produced by the polarized antiproton beam without any free 
parameter. The comparison of model prediction with experimental data is 
illustrated in Fig. 3. The solid lines are model predictions. The experimental 
data for pp-coUisions are shown by the closed circles, while data for pp colli- 
sion are indicated by the open circles. Relation is consistent with data 
in general (Fig. 3) but there are the drastic discrepancies in some cases which 
are not yet understood. For example, Fig.3a shows, that the experimental 
data violate the rule ()29|) . while model gives a fairly well description of the 
SSA for 7r+ production. In Fig. 3b the model (jH)) gives a better description 
of data than parametrization (0). The data in Fig.3c are consistent with 
expectation and the CMSM gives a good description of SSA. 

The Fig. 3d, 3e, and 3f show the pt dependencies of the SSA for 7r^'° 
inclusively produced by protons and antiprotons. In the explicit form such 
dependencies were presented in the original papers only for reactions p(t)p 
n'^X [13] (Fig.3d and Fig.3e, open circles). Other figures were restored by us 
using the tables of data in the original papers ^Hl, and [121. function 
of Pt all SSA are consistent with the charge conjugation rule (compare the 
closed and open circles in Fig. 3d, 3e, and 3f). 

The solid lines are predictions of the CMSM. These lines were calculated 
by using the relation xp ^ a + bpx deduced from the tables of data in 
the original papers. Constants a and b were defined by fit for each pair of 
SSA as they are presented in figures. The consistencies between data and 
predictions are reasonable, though in Fig.3d and Fig.3e the discrepancies are 
seen at low px between model prediction and the experimental data. At such 
low Pt some non-perturbative interaction ( in particular, instanton induced) 
may affect our predictions by changing the weighting factor w(x) and/or by 
flipping quark helicity. 

Concluding this section we note that in general the CMSM gave the qual- 
itative description of most of the SSA effects in a whole measured kinematic 
domain. Nevertheless it fails in the interpretations of the one feature of SSA, 
namely, the evidence of A^- ^ in the interval < xp < 0.4. This is illus- 
trated by the An{p{^)p +X) (see Fig.3b) and AAr(p(|')p ^ 7r"'"+X) (see 
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Figure 3: Test of the charge conjugation rules (fTTj) and for SSA in 
inclusive pion productions in (closed circles) and p(t)]?-(open circles) 

collisions at 200 GeV/c. SSA are given in function oi xp in Figures a), b), 
c) and in function of px in Figures d), e), and f). Solid lines are the CMSM 
predictions for the initial quark polarizations (jSj) and the dashed lines are for 
conditions (P). 
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Fig. 3d). This event might be related to some threshold effect like resonance 
production. 

Conclusion 

We propose a Chromo-Magnetic String Model producing the simple ana- 
lytic expressions for the single spin asymmetries of inclusive charged and 
neutral pion productions. It describes all known experimental data in the 
central and beam fragmentation regions. The following experimental facts 
are naturally explained on the basis of the parton model, conservation of 
the charge and the invariance of the parton interaction with respect to the 
charge conjugation: 74jv(p(|) + p ^ n"^ + X) ^ An{p{^) + p ^ + X) 
and An{p{]) + p ^ tt^ + X) ^ An{p{]) + p ^ tv^ + X). The qualitative 
descriptions of the pt, xf and energy dependencies of the single spin asym- 
metry seem reasonable. The CMSM should be applied beyond the domain 
of resonance production and out of the threshold phenomena. 
To improve the model one has to use the quark and gluon distributions 
{V{x)) given by the global parton analysis and the fragmentation functions 
{D[z)) measured in e~e"''-annihilation together with the NLO expressions for 
the hard parton-parton scattering cross sections. This will provide a more 
precise weighting factor w{x) ((Tj). 

Next (and more important) is to use a more realistic formula for the quark po- 
larization Pq{x). Unfortunately we can not use the spin dependent structure 
function g2{x) measured in the deep inelastic scattering (DIS) directly. This 
function has no parton interpretation and for 100% polarized on mass-shell 
quark g2{x)=0. Therefore, the value of g2ix) measured in DIS is expected 
to be smaller than Pq{x). On the other hand it looks reasonable to choose 
Pq{x) = c ■ g2{x) with such a coefficient 'c' that (according to "spin crisis") 
the whole polarization carried by valence quark will be about 1/4 of the in- 
coming proton polarization. In this way, thanks to the x behavior of g2{x), 
we hope to reproduce the behavior of SSA observed in the data at xp < 0.3. 
We expect that the coming soon data from the polarized RHIC will test our 
model more precisely. 

We would like to express our thanks to L. Bland, Yu.V. Kharlov, V.V. 
Mochalov and D.V. Sidorov for help and useful discussions. 
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